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ABSTRACT 



(D 

The prompt (i<,0.16 days) light curve and initial 9-th magnitude optical flash 

^ ■ from GRB 990123 can be attributed to a reverse external shock, or possibly to 

f^ . internal shocks. We discuss the time decay laws and spectral slopes expected under 

^r^ , various dynamical regimes, and discuss the constraints imposed on the model by 

C^l \ the observations, arguing that they provide strongly suggestive evidence for features 

Qv^ \ beyond those in the simple standard model. The longer term afterglow behavior is 

^* ' discussed in the context of the forward shock, and it is argued that, if the steepening 

i-C ■ after three days is due to a jet geometry, this is likely to be due to jet-edge effects, 

I ' rather than sideways expansion. 

o 

J^ ' Subject headings: gamma-rays: bursts — shocks — optical radiation - cosmology: 

miscellaneous 

't 

1. Introduction 

The observations of GRB 990123 (Akerlof, etal, 1999b; Kouveliotou, etal 1999; Kulkarni et. 
al 1999) not only pose constraints on the amount of gamma-ray beaming needed from a stellar 
mass progenitor in the absence of lensing, but also provide an interesting test of the canonical 
fireball shock afterglow model. A simultaneous optical flash of 9-th magnitude from a burst at 
cosmological redshifts was discussed more than two years ago by Meszaros & Rees 1997 (models 
a2, a3), from the reverse shock that acompanies the blast wave. This early optical flash is expected 
to start at the time of the gamma-ray trigger, and to decay faster than the better known radiation 
from the forward blast wave, which starts out weaker but dominates the longer duration optical 
afterglow. A similar prediction based on the reverse shock was made in the light of more recent 
studies by Sari & Piran (1999a), and more specifically discussed by them in the context of the 
observations of GRB 990123 (Sari & Piran 1999b). 



A different origin for a simultaneous optical flash is possible from internal shocks (Meszaros & 
Rees, 1997, e.g. model b2 of that paper). Internal shock optical flashes are of additional interest 
because, as pointed out by Fenimore et al (1999), the gamma-ray light curve of GRB 990123 (as 
well as those of several other bursts) appears to be incompatible with the gamma-rays coming 
from a single external shock, since the gamma-ray pulses in the second half of the burst are not 
appreciably longer than in the first half. It is, of course, possible that the gamma-rays arise in 
internal shocks, which are thought to be exempt from such problems, while the optical afterglows 
may arise from the external and the associated reverse shocks. However, in the light of the need for 
internal shocks, it is interesting to investigate the implications of the early afterglow observations 
at various wavelengths including optical, for both external and internal shocks. 

We also discuss the longer term behavior of the afterglow, and the causes for the optical light 
curve flattening after 0.16 days. The likely interpretation of the light curve after this time is that 
it is due to the forward shock or blast wave. The discrepancy pointed out between the observed 
time decay slope and the spectral index within the context of the simple standard model can be 
resolved by invoking the simplest realistic extensions to this model, and we discuss several specific 
possibilities (§4). We also indicate that the steepening of the light curve to oc t~^'^ after about 
three days, if due to a jet geometry, is likelier to be due to the effects of beginning to see the edge 
of the jet. This effect occurs before, and its effects fit the steepening better, than the alternative 
sideways expansion interpretation (§5). 



2. Optical flash from reverse external shocks 

The reverse shock acompanying the forward blast wave gives the right magnitude prompt 
optical flash with reasonable energy requirements of no more than a few 10^^ erg isotropic 
(Meszaros & Rees 1997). The time decay constants calculated in that paper for models a2, a3 
were affected by an error, which we correct here; we consider also a more generic prescriptin for 
the dynamics and the magnetic field behavior. 

For a general evolution of the bulk Lorentz factor F oc r~^ with radius, the radius and 
observer time t are related through r '-^ ctF^, or 

F oc r-5 oc t-5/(i+29) , r oc t^^^^^+^a^ . (1) 

In the usual Blandford-McKee (1976) impulsive solutions, the "adiabatic" case is 
g = 3/2, F oc r~^'^, r oc t^'^; the "radiative" case is g = 3, F oc r^^, r oc t^'"^ and in 
the similarity limit g = 7/2, F oc t^"^'^, r oc t^'^. More general values of g occur if the injection 
is non- uniform (Rees & Meszaros 1998), anisotropic or the external medium is inhomogeneous 
(Meszaros , Rees & Wijers 1998). Strictly speaking F is the bulk Lorentz factor of the forward 
shocked material, and may be used also for the contact discontinuity. It is only a very rough 
approximation for the Lorentz factor of the reverse shock. Assuming the latter approximation is 
valid, the comoving width, volume, and particle density of the ejecta, after it has been traversed 



by the reverse shock, evolve with 



AR ~ r/T oc r^+s , V' ~ n'^/ oc r^AR oc r^^^. (2) 

We consider, in the reverse shocked gas, two possibihties for the comoving magnetic field evolution: 
one is flux- freezing, B' oc y'~^/^, and the other is that the comoving field in the reverse shocked 
gas remains in pressure quilibrium with the forward shocked gas, B' oc T. For flux-freezing 
(pressure eqiulibrium) we have then 

B' oc r-^^+^s^/^ , (or oc r~S) . (3) 

The energy density e' and the electron random Lorentz factor 7 in the reverse shocked gas follow 
from e' oc V'~^/^ oc r-(i2+43)/3^ ^^j. ^c T"^ oc r''^9) and 7 oc e' /n'^^ oc r'^^^+a)/^ (oc r^-^). The 
synchrotron peak energy in the observer frame is then 

U^ oc TB'^^ oc r-(12+79)/3 ^ t-{12+79)/{3+69) 

( ocr6-49oct(6-49)/(l+2c,)) ^ (4) 

in these two magnetic field cases. Considering for simplicity the case where the electron cooling 
time is long compared to the dynamic expansion time, the comoving synchrotron intensity at the 
peak frequency is I^^ oc n'^jB'AR oc r^^^'^^'^^''^ (oc r^(^+^^), and the observer-frame flux is 

F,„ oc t^r^'ll^ oc r-(6+59)/3 oc i-(6+59)/(3+69) 

( OCr-29oct-29/{l+29)) _ (5) 

For a photon energy spectral index /? [F^ oc v^) the spectral flux at a given frequency (e.g. optical) 
expected from the reverse shocked gas is then 

^ ^-[6-12/3+5(5-7/3)1/3 ^ ^-[6-12p+g(5-7P)]/{3+6g) 
( 0^^-[23(l-2/3)+6/3] o^^-[2g(l-2/3)+6/3]/{l+29)) ^ (g) 

where /3 = 1/3 below i/m and /? = —{p— l)/2 above i^m in synchrotron radiation. We have assumed 
here that the cooling frequency is above the peak frequency, the conditions for the latter appearing 
to be satisfied in GRB 990123, as pointed out by Sari & Piran 1999b. Under the flux-freezing 
field behavior, for an adiabatic case g = 3/2 and an electron index p = 2.5, the photon spectral 
index above z^m is /3 = —3/4 and we have F^ oc t~8i/32 ^ ^-2-5^ whereas for p = 2 one would have 
/? = —1/2 and F^ oc t~33/i6 ^ ^-2 (^j^j^g fg^. ^j^g similarity case g = 7/2, p = 5/2 we get, as Sari &: 
Piran 1999b, F^ oc t~^^^'^^^ ~ t~^'^). Under the pressure equilibrium field behavior, for g = 3/2, 
and arbitrary p, we have F^ oc t~^'^ (this is a degenerate case where u^ is constant), whereas for 
g = 7/2, p = 5/2 we have F^ oc t"^^/^ ^ t^^-*^. The ROTSE observations (Akerlof etal, 1999) 
give an approximate dependence oc t^^ for about 600 seconds), which is in rough agreement with 
the flux-freezing value in either the adiabatic p = 2 or similarity p = 5/2 cases, which cannot be 



distinguished without spectral information during the ROTSE observations. A different analysis 
of the earliest data on GRB 990123 (Fruchter et a/.1999) gives the ROTSE slope as oc t~^-^ , which 
is close to the pressure equilibrium similarity solution for p = 5/2. 

However, note that a decay oc t^^ (or oc t^^'®) can also be obtained in more generic situations 
than the above, and in particular the spectral slope need not be the only constraint on the decay 
index, except in the simplest, homogeneous external medium, or single initial T (impulsive) model. 
For an impulsive injection in an inhomogeneous external medium, e.g. p^xt of r^ one expects in 
the adiabatic limit F oc r^^ with 5 = (3 — d)/2 (Meszaros , Rees &: Wijers 1998). For a more 
realistic non-uniform injection situation, one expects a range of initial F, and in particular for a 
power-law distribution of F in the ejecta (Rees & Meszaros , 1998) where the mass injection rate 
M(> F) oc F~*, the effective resulting F oc r~^ with g = 3/(2 -|- s). The reverse shock relation 
(6) between photon spectral slope and dynamics could therefore, in the simplest "post-standard" 
extensions of the model, depend e.g. on the external density profile, or the mass-energy injection 
dependence. 



3. Optical flashes from internal shocks 

Internal shocks, just as external shocks, should have a low energy portion or tail of the 
spectrum which extends into the optical, starting at the GRB trigger (except possibly in rare 
cases where the self-absorption frequency extends to the optical). For the same total energy, the 
flux from the internal shock optical flash is about 2 orders of magnitude weaker than from the 
reverse external shock (Meszaros & Rees 1997); however, for beaming factors ~ 10~^ as widely 
considered, they could lead to flashes as bright as 9-th magnitude at z --^ 1. The simple "standard 
internal shock model" (Rees & Meszaros , 1994) considers a wind of duration t^ ~ tiy^^^-gt with 
< F >'^ constant and fluctuations AF ~ F over times t^ < t^, producing shocks at a distance 
rgh ~ ct^, < F >^ for t<tw 1) The most straighforward GRB internal shock (e.g. model bl of 
Meszaros & Rees 1997) assumes that the gamma-ray break observed around ~ 100 keV is due to 
the synchrotron peak energy. The magnetic field required in the internal shocks is not far below 
equipartition, and this leads to a very short radiative cooling time compared to the expansion 
time. If the shocks arise from discrete, sharp-edged shells, there would be a very sudden drop 
of the light curve at all energies (including optical) after internal shocks stop. However, a more 
realistic situation probably involves smoothly modulated shells, where the shocks weaken gradually 
after reaching a peak strength. A variety of plausible distributions would then lead to a power law 
decay of the optical light. 

2) An alternative internal shock model, which would lead to long decay even in the absence of 
a gradual fading of the shocks, is obtained in the case where the synchrotron break is at optical 
energies, and the gamma-ray break is produced by IC-scattered synchrotron photons (model b2 of 
Meszaros & Rees 1997, Papathanassiou and Meszaros , 1996). Between the time at which internal 
shocks (or gamma-rays) stop and the time when deceleration by the external medium starts 



{tdec ^ 500(£'54.6/(^/47r)ne2:t)^'^r3QQ sec), the average bulk Lorentz factor remains approximately 
unchanged, and the dynamics are described by F ex i^, r oc t. Hence, the comoving width, volume, 
and particle density of the ejecta evolve as 

Ai? ~ r/r oc t , y ~ n'^j ~ r^Ai? oc t^ (7) 

and the comoving energy density and electron random Lorentz factor will be 

e' ~ n'7 oc V'"^/^ oc t'^ , 7 ~ e' /n' oc r^ . (8) 

If the magnetic field is not dynamically dominant, two cases are 

a) B' oc V'^^'^ oc t^^, if the field is random, and 

h) B' oc B (X r^^ oc t^^, if the field is mainly transverse (e.g. inefficient reconnection). Cases a (b) 

lead to an observer-frame synchrotron peak flux 

Um oc T.B'.-f^ oc t-"^ (oc r^) (9) 



(The first dependence assumes random fields , and the second dependence is for transverse fields). 
The comoving intensity is then /^^ oc n'^jB'AR oc t~^ (oc t~^) and 

F,^ oc t^r^/;,^ oc t-2 (o,t-i)^ (10) 

while the flux at a fixed frequency is 

where one can verify that at the internal shock radius and above, the optical electrons are in the 
adiabatic regime, for a wind equipartition parameter es = 10~^ and 7 -^ 300. Thus for p = 2, 
F^(xt~^ (oct-5/2) 

For a magnetically dominated outflow, an interesting third case is 
c) B' (X B (X t~^, as in (b) above, but now one might expect the comoving volume and particle 
density to be dominated by the field evolution (via B' oc t~^ oc V'^'^''^), so 

y oc n';/ oc t3/2 ^ AR^V'/r^ocr^/^ (12) 

e' ~ n'7 oc V'~^/^ oc t-2 ^ ^ _ ^'/n' oc r^^^, (13) 

Um OC T.B'.-f^ oc t-2 (14) 
the comoving intensity is I^^ oc n'^jB'AR oc t^^ and 

F,„ cc t'T'll^ cc t-i (15) 

F. - F.^v^r'^'^ - t-^. (16) 

Adiabatic conditions prevail for e^ -^ 1 and, for example, F>300, 7 ^ 100. For p = 5/2 this gives 
Fjy oc t^^i"^, and for p = 2 it gives F^ oc t^"^. 



For the simple internal shock model where the average T is constant, model (a) is too steep, 
but model (b) and (c) could fit the observations, with a flat enough electron power law (p ~ 2). 
However, in a realistic model the average F of the wind producing internal shocks could vary, and 
even for p > 2, a decay t~^ or t~^'^ could be the result, e.g., of an average F which increases in 
time. 



4. The external blast wave 

The standard long-term (as opposed to short-term) afterglow is attributed to the external blast 
wave (forward shock) evolution, whose bulk Lorentz factor is described by the same equations (1), 
and the comoving width and volume of the forward shocked gas by equations (2). The comoving 
density of the forward shocked gas is, however, n' ~ riextF and the comoving field is assumed to be 
some fraction of the equipartition value, B' oc e^ F, while the shocked electron random Lorenzt 
factor is 7 -^ Se^- As in Meszaros & Rees 1997 (model al), Vm of r~^^ oc t~^5/(i+29)^ £qj. adiabatic 
electrons/' oc n'B'AR (x r^~^, F^, (xt^T^I' oc r'^''^^ oc t'-^"'^^^ / '^^^'^^\ andioi a, STieciiwuKx v^ 

Urn ' '^"^ '^m ' ■'- 

we have 

p^ ^ ^3-25(1-2/3) ^ ^[3-29(l-2/3)]/{l+2a)_ (^^7^) 

Thus in the (impulsive) "standard model" a simple relation is expected between the time decay 
index a and spectral slope (3, namely a = [3 — 2g{l — 2/3)]/(l -|- 2g). E.g. for the impulsive 
adiabatic case g = 3/2 and one would expects a = (3/2)/3 with /? = 1/3, — (p — l)/2 for 
synchrotron, or their equivalent for the radiative ^f = 3 or similarity g = 7/2 values. Such a simple 
one-parameter relation between a and (3 does not appear to hold for the second stage of GRB 
990123, where Fopt oc t~^'^ (Kulkarni, et al., 1999). This, in our view, is a strong indication that 
"post-standard" features are present, e.g. a non-homogeneous external medium or an anisotropic 
outflow (Meszaros , Rees &: Wijers, 1998), or non-uniform injection (Rees & Meszaros , 1998). 
From an observed time decay Fopt ^ i~" and an observed spectral slope F,y oc v^ one can then 
work backwards to get the effective value of g implied by the equations (6) above, and this in turn 
implies a value oi d = 3 — 2g or s = {3/g) — 2, thus providing information about the external 
medium or the injection mechanism. While non-unique, such examples indicate that the time 
decay index is likely to depend on parameters other than the spectral slope. 



5. A Jet Geometry? Edge vs. Expansion Effects 

The temporal decay index after three days steepens from about t^^'^ to t^^'^ (Kulkarni, 
et aZ.199; Fruchter, et a/. 1999), and this steepening can be attributed to a collimated outflow, 
whose effects become noticeable after F drops sufficiently. If one assumes (cf Kulkarni, et a/. 1999) 
that the steepening is caused by sideways expansion of the decelerating jet (Rhoads, 1997) one 
expects a large steepening from t~^'^ to i~^, where p is electron index, typically p = 2.5, so the 



change expected would be more than one power of t. However, the edge of the jet begins to be 
seen when T drops below the inverse jet opening angle l/Oj. This occurs well before sideways 
expansion starts (Panaitescu & Meszaros , 1998): the latter is unimportant until the expansion 
is almost non-relativistic. So long as F > ^~ , the emission we receive is the same as from a 
spherically-symmetric source, and the effective transverse area is A ~ (rii/r)^ oc t^F^; on the 
other hand, when F < ^^ the dependence is ^ ~ (^ll^i)^ °*- ^^r^. Note that even after the edge 
of the jet becomes visible, the outflow is still essentially radial and F continues to decay as a 
power law as before (until sideways expansion sets in). These two additional powers of F (e.g. in 
adiabatic expansion, F oc r"^'^ oc t~^'^) imply a steepening by t~^'^, which matches quite well the 
observed steepening by a power of about t^^'^. There are of course several other plausible causes 
for steepening, as discussed by Kulkarni et a/.(1999), Meszaros , Rees & Wijers 1998, Rhoads 1997 
and others. However, the agreement between the observed change in the decay slope and that 
expected from seeing the edge could well be a signature of the detection of a jet. 



6. Discussion 

The conclusions that may be drawn from the above are that bright optical flashes, starting 
at the time of the gamma-ray trigger and extending into the early afterglow, are a robust 
prediction of the simplest afterglow models, as discussed in Meszaros & Rees, 1997, Sari &; 
Piran, 1999a, 1999b. A 9-th magnitude optical flash from a GRB at redshift z ~ 1 can arise 
from a reverse external shock with total energy £'53(i7/47r)^^ ~ 1, or from internal shocks with 
£;53(ri/10"247r)~i(^/102s)"^ ~ 1 (Meszaros & Rees 1997), where Q is the solid angle into which 
the radiation is collimated, and t^ is the GRB duration. As known from logA^ — logP fits (e.g. 
Krumholz etal, 1999), in essentially all cosmological models the luminosity function must be broad, 
with a range of E which can span upwards of 2-3 orders of magnitude. This is compatible with 
the fact that GRB 990123 is in the top 1% of the BATSE brightness distribution (Kouveliotou 
etal, 1999), and also with the previous non-detection of similarly bright optical flashes. 

A steeper time-decay law is expected for the early optical flash from reverse or internal 
shocks, compared to that expected from the forward shock. The decays calculated in Meszaros & 
Rees 1997 are too steep, compared to the observed t~^ behavior (Akerlof etal 1999), indicating the 
need for investigating different assumptions for the magnetic fleld behavior. In the present paper 
we find that an early time-decay oc t^^ is naturally explained in the standard (single-F) afterglow 
model either by reverse shocks, where the field is in pressure equilibrium with the forward shock 
or is frozen-in, or by internal shocks in a magnetically dominated outfiow. However, such decays 
may also be easily achieved in more realistic afterglow models without severe restrictions on the 
electron index, e.g. for bursts in an inhomogeneous external medium, or characterized by a power 
law range of F. 

From the fact that the inferred isotropic equivalent gamma-ray energy is 4 x 10^ ergs 
(Kouveliotou etal, 1999), one infers the need for a collimation of the gamma-rays by at least 10~^, 



if not 10^^, for any stellar-mass source. For an external reverse shock origin, from the energy 
alone the optical radiation need not be collimated at all; however, since in the observer frame the 
Lorentz factor of the reverse shock is initially close to that of the blast wave, the initial optical and 
gamma-ray collimation should be approximately the same (except, as may be the case, if the jet 
has an anisotropic T). For an internal shock origin of the optical flash, there are stronger grounds 
for expecting the beaming of the optical to be the same initially as for the gamma-rays. 

The flattening of the light curve after 0.16 days to oc t^^'^ is well explained by the optical 
light expected from the forward shock, conforming to the standard interpretation of the long-term 
behavior of afterglows. The discrepancy pointed out Kulkarni et a/. 1999 between the standard 
model prediction and the observed relation between the time decay index and the spectral slope 
is, as pointed out in §4, most likely to be an indication of departures from the simple standard 
model, e.g. an inhomogeneous external medium, non-uniform injection, or anisotropies. The 
renewed steepening of the optical light curve after about three days to oc t^^'^ (Kulkarni, et 
aL1999, Fruchter, et a/.1999) is, as argued in §5, likely to be due to seeing the edge of a jet, which 
gives a better fit for the magnitude of the change of slope than expected from sideways expansion, 
and also should occur well before the latter begins. 

There are two main arguments that militate in favor of the optical flash observed in GRB 
990123 being due to a reverse external shock. One is that the optical light curve does not show a 
good correlation to the gamma-ray light curve (Sari & Piran ,1999b), and the model fits generally 
well. The other is that the fiash from internal shocks is weaker than from external reverse shocks, 
as discussed here. On the other hand, the first ROTSE observations of GRB 990123 (Akerlof 
et al, 1999a) started 20 seconds after the GRB trigger, and do not appear to have sampled the 
optical light curve densely enough to establish the degree of correlation with good significance 
(both curves are compared in Fenimore et al., 1999). While the likelihood of detecting an external 
reverse shock optical fiash appears to be higher at this stage, one needs a faster triggering and a 
better sampling of the optical data in order to discriminate between reverse and internal shocks. 
An interesting possibility is that, since the external reverse shock starts somewhat later than the 
internal shocks, one might initially see a weaker optical fiash from the internal shocks, which are 
overtaken by a stronger reverse external shock radiation after tens of seconds, until the forward 
shock optical afterglow takes over after 300-1000 s. Especially for the more frequent weaker bursts, 
the prospect of investigating such features underlines the need for dedicated GRB afterglow 
missions such as HETE2 and Swift. 
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We are grateful to A. Panaitescu for useful insights. 
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